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ABSTRACT 

We have used the IRAM Plateau de Bure Interferometer and the Expanded Very Large Array 
to obtain a high resolution map of the CO(6-5) and CO(l-O) emission in the lensed, star-forming 
galaxy SMM J2135— 0102 at z = 2.32. The kinematics of the gas are well described by a model of a 
rotationally-supported disk with an inclination-corrected rotation speed, v Iot = 320±25kms _1 , a ratio 
of rotational- to dispersion- support of vja = 3.5 ± 0.2 and a dynamical mass of (6.0 ± 0.5) x 10 10 M Q 
within a radius of 2.5 kpc. The disk has a Toomre parameter, Q = 0.50 ± 0.15, suggesting the 
gas will rapidly fragment into massive clumps on scales of Lj ~ 400 pc. We identify star-forming 
regions on these scales and show that they are ~ 10 x denser than those in quiescent environments 
in local galaxies, and significantly offset from the local molecular cloud scaling relations (Larson's 
relations). The large offset compared to local molecular cloud lincwidth-size scaling relations imply 
that supersonic turbulence should remain dominant on scales ^100x smaller than in the kinematically 
quiescent ISM of the Milky Way, while the molecular gas in SMM J2135 is expected to be ^50x denser 
than that in the Milky Way on all scales. This is most likely due to the high external hydrostatic 
pressure we measure for the interstellar medium (ISM), Ptot/^B ~ (2± 1) x 10 7 K cm~ 3 . In such highly 
turbulent ISM, the subsonic regions of gravitational collapse (and star-formation) will be characterised 
by much higher critical densities, n cr i t >=10 8 cm -3 , a factor >1000x more than the quiescent ISM of 
the Milky Way. 

Subject headings: galaxies: starburst, evolution, high-redshift, gas, star-formation; galaxies: individ- 
ual: SMMJ2135 0102 



1. INTRODUCTION 

Observations of the molecular gas in galaxies provides 
important insights into the physics of star formation. 
Since all stars form out of molecular gas - either in self- 
gravitating giant molecular clouds (GMCs) in disks such 
as the Milky Way or distributed in continuous gaseous 
disks in ultra-luminous in frared galaxies (e.g. ULIRGs; 
iDownes fe Solomon! I1998T ) - studying its thermal and 
kinematic state is a prerequisite for understanding star 
formation through-out the Universe. Such observations 
of the interstellar medium (ISM) , especially in the most 
extreme conditions like those in the high star-formation 
density environments of local ULIRGs, can provide pow- 
erful tests of star-formation theories and the initial con- 
ditions of the stellar initial mass function (IMF). In 
the most energetic systems, strong star-formation feed- 
back, high cosmic-ray energy densities and supernovae- 
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driven shock heating of very turbulent molecular gas al- 
lows us to explore uncharte d parameter space in cur- 
rent star-formation schem es (|Krumholz fc McKed 120051 
iPapadopoulos et all 120 111) . 

Given that most of the star formation in the most mas- 
sive galaxies is claimed to have occurred at z ~ 1-3, an 
era when the gas accretion rate and molecular gas frac- 
tions o f galaxies appear to be substantially hi gher than 
today (|Tacconi et aijl201Ct iDaddi et all [2010). examin- 
ing the physical and dynamical states of the molecular 
gas in a star-forming system during that era acquires 
special importance. Evidence is accumulating that the 
dominant mode of star formation may be very different 
in early systems than that found in most local d isks (e.g. 
iBournaud fc Elmegreenll2009l I Jones et al.ll2010h . Rather 
than forming stars within GMCs that condense out of a 
stable disk, star formation could be triggered by fragmen- 
tation of dynamically unstable gas-rich disks. The result- 
ing regions of high-pressure ISM may be conducive to the 
for mation of massive g lobula r clusters, as first postulated 
by lHarris fc Pudrita ()1994[ ). a model which has since 
been developed extensively (e.g. lElmegreen fc Efremqvj 
[19971 ISchweizer etldl [19961 IKravtsov fc Gnedinl 12005). 
Such a star-formation mode likely occurs locally, albeit 
on more compact scales, in the most extre me merger- 
driven starbursts (jDownes fc Solomonlll998H . A spatial 
resolution of < 100 pc is necessary to resolve the ex- 
pected size of gravitationally unstable regions in gas-rich, 
star-forming disks, and to probe the masses, sizes and 
hence densities of the molecular ISM on these scales. 
For example, the estimated size of the two gas-rich, 
star-forming disks of Arp 220 - an archetypal ULIRG 
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in the local Universe - is ~ 100 pc, as only recently re- 
veale d by sub-millimetre interferometry (Sakamoto et al.l 
2008). Achieving this spatial resolution within starburst 
galaxies at z ~ 2 (~ 200x more distant than Arp220), 
with sufficient sensitivity to measure the properties of 
the molecular ISM, is one of the key science drivers for 
the Atacama Large Millimetre Array (ALMA). However, 
prior to full science operation of ALMA, progress can 
still be made by exploiting rare examples of high-redshift 
starburst galaxies which have been fortuitously strongly 
gravitationally lensed by foreground clusters, boosting 
both the apparent size and flux of the galaxy and so al- 
lowing imaging spectroscopy at high spatial resolution 
with current millimetre arrays. 

The recent discovery of a gas-rich starburst galaxy 
at z = 2 3, SMMJ213511.60- 010252.0 (hereafter 
SMMJ2135; iSwinbank et al J 120 lOf ) . provides just such 
an opportunity. The apparent brightness of SMM J2135 
is due to a factor ~35x magnification by the fore- 
ground, mass i ve g alaxy cluster, MACS J2135— 0102 
(jEbeling et al.l 120011) . This has allowed very detailed 
studies of the distribution and intensity of star formation 
and the chemical make-up of the ISM in this example of 
the cos mologically important high-redshift ULI RG popu- 
lation ([Ivison et al.ll2010l : iDanielson et al.ll2011l ). In par- 
ticular, high-resolution interferometry with the Smith- 
sonian Sub-millimeter Array (SMA) resolved rest-frame 
260-/im continuum emission in two mirror-image sets of 
four clumps straddling the critical line. The clumps have 
physical (source-plane) sizes of ~ 100-200 pc and are dis- 
tributed over a region ~ 2 kpc within the background 
galaxy. 

If we can similarly resolve the ISM on 100-200 pc 
scales, this would allows - for the first time - a direct 
comparison of the scaling relations found i n local molec - 
ular clouds: the so-called Larson relations (jLarso n 1981). 
These relate the CO velocity line width of clouds (a) with 
their physical extent, R (a oc R p , p ~ 1/2) and mean 
molecular gas density with size ((^(Hb)) oc R~ k , k ~ 1). 
These relations have been found to be valid through- 
out the Milky Way GMCs as well as in those extra- 
galactic environments whe re high-resolution s tudies of 
GMCs have been possible (|Bolatto et al.l 120081) . In par- 
ticular, the universality of the a-R relation within star- 
forming GMCs points towards large-scale driving mecha- 
nisms, rather than small-scale turbulent energy injection 
from within molecular clouds, as the origin of interstel- 
lar turbulence. The l atter is a driver of th e properties of 
star- forming clouds (jHever fc Bruntl 120041 ) and may de- 
fine the star-formation efficiency and the result i ng stel- 
lar Initial Mass Function (IMF) (iKlessenl l2004t iLarsonl 
120051 Uappsen et "all 120051: iKlessen et al.ll2007f ). In this 
context, Larson's relations provide a unique probe of the 
dynamical state of the turbulent molecular gas in extra- 
galactic star-forming systems, which are well-calibrated 
locally and are complementary to molecular line ratios 
that constrain the density and temperature (i.e. the ther- 
mal state) in these environments. 

In order to understand the physical processes that 
drive the apparently high star-formation rates within 
compact star-forming regions at high redshifts, we 
have obtained high-spatial-resolution observations of the 
molecular gas in SMMJ2135 traced by the CO(6-5) 



emission with the IRAM Plateau de Bure Interferometer 
(PdBI). We use this spatially- resolved CO spectroscopy 
to probe the dynamics of the galaxy, and the luminosities 
and dynamics of the star-forming gas. We independently 
confirm our results using a comparable resolution map of 
the CO(1-0) data taken with the Expanded Very Large 
Array (EVLA). In §2 we detail the observations and data 
analysis. In §3 we discuss the properties of individual 
star-forming regions. In §4 we discuss the wider context 
of our results for star-formation theories, and we present 
our conclusions in §5. Throughout the paper we use a 
ACDM cosmolo gy with H n =72km s~ 1 Mpc" 1 , fi m =0.27 
and n A =l-Vt m ([Spergel et all 120031 ) 



2. OBSERVATIONS 

We obtained Hubble Space Telescope (HST) WF3 
F110W- and F160W- band images (each of exposure time 
1.2 ks) on 2010 August 23 as part of our snapshot pro- 
gram GO# 12166; PI: H. Ebeling). The images com- 
prise three 400 s exposures dithered with a 3" spacing. 
The data were reduced using MULTIDRIZZLE to provide 
an image with 0.12" sampling and good cosmetic prop- 
erties. A visual inspection of the image clearly iden- 
tifies the red, radial image of SMMJ2135 located pre- 
cisely at the position identified by the low resolution 
discovery- and high-resolution follow-up- observations at 
a : 21:35:11.6 5 : -01:02:52.0 (J2000) with a spatial ex- 
tent of ~7.5". A search of the WF3 image also reconfirms 
the counter-image of SMM J2135-0102 at a : 21:35:15.6 
S : -01:03:13.12 (J2000). In Fig. Q] we show the observed 
VJJ J"-band colour colo ur image around the cluster core 
(see iSmail et al.l 120071 for a description of the HST V- 
band observations). We also highlight the position of 
a second z — 2.3 triply-imaged system (labelled arc2- 
a/b/c) located ^37" due South of the cluster core. We 
will return to a discussion of this imaging data, in par- 
ticular in the context of the lens modelling, in § 13.11 

The redshifted CO(6-5) emission from SMM J2135 was 
observed with the six-element IRAM Plateau de Bure 
Interferometer (PdBI) at 207.90 GHz in "A" configura- 
tion. We obtained a total on-source observing time of 
6hr on 2010 January 6 . In this configuration the syn- 
thesised beam for natural weighting is 0.67 x 0.43arcsec, 
at a position angle (PA) of 24.5°. The spectral cor- 
relator was adjusted to give a frequency resolution of 
2.5 MHz across the 980-MHz bandwidth. The overall flux 
density scale was set on MWC349, with observations of 
2134+004 for phase and amplitude calibration. Receiver 
bandpass (RF) calibration was performed on 1749+096. 
The data were calibrated, mapped and analysed in the 
GILDAS software package. In the map, the rms is 1.5 mJy 
per 25-kms -1 channel. Fitting beams to the velocity 
integrated cube, we measure a velocity integrated flux 
of lg+Uykms " 1 , thus recovering > 90 per cent of the 
single dish flux (jDanielson et al.ll2011[ ). 

In addition, the redshifted CO(1-0) line emission was 
observed with the EVLA between 2010 September and 
2011 February, as part of the Open Shared Risk Ob- 
serving (OSRO) period. We programmed the WIDAR 
correlator to return two sub-bands of 64 x 2-MHz dual- 
polarisation channels. To avoid centering the line in the 
noisy channels at the edge of each sub-band, we shifted 
the centroid into the 40 th channel of the lower sub-band, 
and allowed for 10 channels of overlap between the sub- 
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Fig. 1— Hubble Space Telescope (HST) ACS/WF3 V 6 06>/ll0#i40-band image of the massive cluster MACS J2135-0102. We overlay the 
2 =2.32 critical curve from the best-fit lensing model (red solid curve), and also identify the image of SMM J2135-0102, as well as the 
triply-imaged galaxy at z =2.3 (labelled arc2-a/b/c). The crosses show the cluster galaxies which are included in the lens modelling. The 
insets show: top: a zoomed image of SMM 2135-0102, middle: zoomed image of SMM J2135 with the CO(1-0) emission from the integrated 
cube as contours (marked at 3,6,9...fx); and bottom: the counter-image of SMM J2135-0102 located approximately 53" due East of the 
Brightest Cluster Galaxy (BCG). The co-ordinates for the image (and insets) are centered on the BCG at a :21 35 12.12 8 -01 02 58.80 with 
North Up and East Left. 



bands, yielding a total dual-polarisation bandwidth of 
236 MHz, or 2,040 km s _1 of velocity coverage across the 
line. 

As a result of dynamical scheduling, data were ob- 
tained in excellent weather conditions: yielding seven 2- 
or 3-hr blocks, totaling 10 and 5 hr in the hybrid DnC and 
CnB configurations, respectively, although with a similar 
data volume in each case. During the blocks we switched 
every few minutes between the target and the nearby, un- 



usually bright (4-5 Jy), compact calibrator, 2136+004, in 
order to track changes in phase and amplitude, and to 
determine the spectral response (bandpass). 3C48 was 
used to calibrate the flux density scale. 

Typically 16 receivers were operational during a CnB 
block after on-line and manual flagging. These data 
were pro cessed using s t andar d ALVS procedures, as out- 
lined by llvison et al.l (|2011[ ) . Imaging of the DnC and 
CnB data was accomplished using IMAGR, yielding a 
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Fig. 2. — Velocity-integrated image-plane maps of the (a) CO(6-5) and (b) CO(1-0) emission in SMMJ2135 from our observations 
with PdBI and EVLA respectively. In both panels, the contours start at 3 a and are spaced by 3 a thereafter, with negative contours 
at the same levels denoted by the dashed lines. The synthesised beam for each observation is shown, bottom right, in each panel, and 
the panels are centered at a: 21:35:11.638 <5: -01:02:52.39 with North up and East left. The red solid line denotes the z =2.32 critical 
curve from the best-fit lens model. Both CO(6-5) and CO(1-0) are extended on spatial scales of ~6arcsec, with the emission mirrored 
across the critical curve, (c): the image-plane CO (6— 5) emission-line morphology, with the rest-frame 260- (im continuum from our SMA 
observations overlaid as contours (Swinbank et al. 2010). d & e: image-plane position-velocity diagram of the CO(6-5) and CO(l— 0) line 
emission extracted from the major axis at a PA 45° east of north (as highlighted by the dashed curve in panel (b). These clearly shows 
velocity gradients of ~500kms _1 across ~6arcsec. The histograms show the collapsed, one-dimensional spectra. T horizontal dotted lines 
denote the approximate positions of the star-forming regions seen in the rest-frame 260/^m cont inuum emission, whilst the arrows (and 
labels) denote the velocities of the four kinematic components identified by Danielson et al. (2011). (f) shows the lens amplification profile 
across the galaxy image for the major axis cross-section defined by the black dashed line shown in panel (b). This plot shows that the 
south-eastern image of the lensed image is less highly amplified than the north-western image (by an average factor ~1.4), which gives rise 
to the observed flux ratio of the northern/southern images in the sky-plane. 



~ 0.81 x 0.60-arcsec synthesised beam (for ROBUST=5) 
at PA = 84°. In the map, the rms is 0.13 mJy per 25- 
kms -1 channel. In this cube, we measure a velocity inte- 
grated CO(1-0) line flux of of 2.1Jykms~ 1 , thus recov- 
ering > 98 per cent of the peak single-dish flux density 
(jSwinbank et al.ll2010D . 

3. ANALYSIS AND DISCUSSION 

3.1. Lens modelling 

Since a necessary component in interpreting the CO 
kinematics of SMMJ2135-0102 is the reconstruction of 
the source, we first review the lensing model constraints. 
The primary constraints on the strong lensing model 
for MACS J2135-0102 are the redshift and location of 
the three images of SMMJ2135-0102 as well as those 
of another spectroscopically confirmed triply imaged 



galaxy at z = 2.3 located approximately 37" due south 
of the brightest cluster galaxy (labelled arc2-a/b/c in 
Fi g. [B . A lensing mode l for the cluster was presented 
in iSwinbank et all (|2010f) . However, with the spatially 
resolved spectroscopy provided by our new observations, 
we can improve the strong lensing constraints using the 
kinematics of SMMJ2135-0102 to precisely locate the 
z=2.32 critical curve. 

We begin by constructing the velocity-integrated emis- 
sion line maps from the CO(l-O) and CO(6-5) observa- 
tions, and as Fig. [5] shows, which are extended on scales 
of 4-6" (see also Fig. [3]) . We als o overlay the observed 
870-^im continuum emission from ISwinbank et al.l (2010) 
which highlights that the positions of the brightest re- 
gions in the sub-millimetre continuum appear to broadly 
align with the molecular gas emission. To investigate the 
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Fig. 3. — Left: Volume-rendered image plane cube of the CO(6-5) emission. The cube is colour coded by velocity, and volume rendered 
according to the intensity in each pixel. The 'fold' in the velocity field denotes the location of the critical curve where the image is mirrored. 
The SMA 870/jm (rest-frame 260^tm) is projected on to the back face of the cube and we also include a faint grey surface which denotes 
the critical curve from the best fit lens model. Although the velocity integrated map in Fig. [2] appears relatively smooth, this rendered 
image clearly suggests that when viewed in velocity space, the gas has a "clumpy" morphology. The full movie of SMM J2135-0102 can 
be downloaded from http://astro.dur.ac.uk/~ams/PdBImovie/ Middle Position velocity diagram for CO(6-5) extracted from the major 
kinematic axis (PA 45" east ot north, as in fig. [^J . Here, we overlay the model rotation curve on both mirror images, and use this to extract 
a luminosity profile of the CO. Right: The CO luminosity profile extracted from the PV diagram to emphasise the clumpy morphology 
of the CO. So that the structure in the CO can be seen against the noise, the error-bars are generated by converting the CO flux density 
into number of photons detected per channel (accounting for the efficiency of the telescope). We show both the CO(6-5) (black) and the 
CO(l-O) (blue) (a detailed comparison between the cold molecular gas seen through the CO(6-5) and CO(l-O), radio and sub-mm emission 
at high resolution will be presented in a forthcoming paper; R. Ivison 2011 in prep). We also highlight the positions of the SMA clumps 
as solid bars which shows that the high star-formation density regions seen at rest-frame 260/xm coincides with the clumps observed in the 
cold molecular gas. 



gas dynamics of the galaxy, we extract position-velocity 
(PV) diagrams from the major axis of the image-plane 
CO(6-5) and CO(l-O) cubes (at a position angle, PA 
= 45°, east of north as highlighted by the dashed line 
in Fig. [2| and show these in Fig. [3] These PV dia- 
grams show that the galaxy has a velocity gradient of 
~ 500 kms -1 across ~4-6arcsec which is folded (mir- 
rored) about the critical curve. The sharp rise and flat- 
tening of the CO emission in the PV diagram has the fa- 
miliar kinematic signature of a rotating system, and we 
discuss the dynamical properties in § 13.21 The dynamics 
of the galaxy in the image plane allow us to kinemati- 
cally locate the position of the critical curve, and by fold- 
ing and cross correlating the velocity field with itself, we 
kinematic ally locate it 0.3±0.1 // f urther south than first 
chosen by fSwinbank et alj (|2010l ) (where it was identi- 
fied from high resolution imaging alone). We therefore 
update the lensing model accordingly and briefly discuss 
the lens model here. 

The cluster lensing model for M ACS J2135-0102 is 
derived using th e LENSTOOL code (jKneib et all 119961: 
Uullo et all 120071 ) whic h employs a parametric mo del of 
the mass distribution (jRichard et all [20091 [2010h . We 
use a model with a single cluster-scale mass compo- 
nent (dark matter halo), as well as individual galaxy- 
scale mass components centered on each cluster mem- 
ber (selected from their V — I colours and shown 
on Fig. [T|). We parametrise the galaxy-scale mass 
components as truncated pseudo-isothermal elliptical 
mass distributions, each described b y the following pa- 
rameters: {%, y, e,0, r cor e , r cu t, fdispl (|Kassiola fc Kovnerl 
Il993t [Kncib ct al. 1996) which denote the the [x/y] cen- 
ter, ellipticity, position angle, core radius, cut-off radius 
and velocity dispersion respectively. To reduce the num- 
ber of free parameters and match the constraints above, 



the geometrical parameters (x, y, e, 9) describing the clus- 
ter galaxies (including the brightest cluster galaxy) are 
matched to their 2-dimensional light distributions. We 
assume that the cluster galaxies follow a scaling relation 
with constant mass-to-light ratio according to an L* clus- 
ter galaxy, and their re spective r nn re and tw parameters 
are scales accordingly (jRichard et al J 1201 01) . Turning to 
the parameters describing the cluster-scale mass, we fix 
r cut at 1 Mpc which is motivated by its very small influ- 
ence on the location of the critical line at large values. 
This leaves [x,y] e, 9, r core and Vdi sp as variables. 

The primary constraints in defining an acceptable 
strong lens model for cluster lensed are that it reproduces 
the locations and redshift of the sp ectroscopically con - 
firmed triple- image lens systems (e.g. lSmith et al.1 12009'). 
and in the case of MACS J2135-0102 we also demand 
that it precisely predicts the location of the critical curve. 
The best-fit model (Table 1) obtains a fit with an r.m.s. 
of 0.15" between the predicted and observed positions 
of the multiple images. Since the lenstool code in- 
corporates a Markov Chain Monte Carlo (MCMC) sam- 
pler, we also generate 1000 error models that lie within 
the 68% confidence intervals around the best fit pa- 
rameters (after marginalising over the other parame- 
ters in each case), and construct the amplification maps 
in each case. In the best fit model, the cluster-scale 
component is centered South East of of the BCG (Ta- 
ble 1), and the enclosed mass within an aperture of 
250 kpc is M=3.3±0.3xl0 14 M Q with an Einstein radius 
of 6» B =34.5±2.0" for z=2.32 (Richard et al. 2011 in 
prep). In Fig.[T$£[2]we show the critical curve on the HST 
and CO imaging, as well as amplification as a function of 
position across the long axis of the arc (and the associ- 
ated lcr error from the family of models which adequately 
describe the lensing mass). The updated model sug- 
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gests a luminosity weighted amplification of the system of 
/z=37.5±4.5 which higher, but within la o f the p revious 
estimate of 32.5±4.5 from iSwinbank et al.l l|2010| ). 

3.2. Gas distribution and dynamical properties 

To construct a two-dimensional map of the velocity 
field from the CO emission in the source plane we use the 
best-fit lensing model to calculate the mapping between 
image- and source- plane, and ray-trace each pixel to re- 
construct the source-plane datacube. We then fit the 
CO (6-5) emission at each pixel with a Gaussian profile 
using a \ 2 minimisation procedure, accepting a fit only 
if the emission is detected at > 5a. Where a fit is made, 
we measure the central velocity, line width (a) and line 
flux. The resulting source-plane velocity field of the gas 
is shown in Fig.Uand resembles a rotating system with a 
peak-to-peak velocity of 240 ± 25 km s _1 within a radius 
of ~ 3kpc. In Fig. |4j we also show the extracted one- 
dimensional rotation curve of the CO(6-5) emission as 
well as the line-of-sight velocity dispersion. 

In order to estimate the disk inclination and true rota- 
tional velocity, we model the velocity field with a rotating 
disk. We use an arctan function to describe the shape 
of the rotation curve such that v(r) = v c &rctan(r /r t ), 
where v c is the asymptotic rotational velocity and rt 
is th e effective radiu s at which the rotation curve turns 
over ()Cour tcau 1997|). We construct the two-dimensional 
kinematic model for the galaxy using six free parameters: 
v c , r t , x, y (the central position), PA and disk inclina- 
tion, i, and use a genetic algorithm with 10 5 random 
initial values to search for a best fit. We demand > 30 
generations are performed before testing for convergence 
to a solution. The best-fit kinematic model is overlaid 
on Fig. H and has r t = 0.4 ± 0.1 kpc, i = 60 ± 8° and 
v t = 275 ± 20kms~ 1 , providing an estimate of the dy- 
namical mass of Mdyn = (6.0 ± 0.5) x 10 10 M for an 
inclination of 60 ± 8° . We note that the total gas mass 
from a large velocity gradient (LVG) analysis, Af gas = 
4 x 10 10 M (jPanielson et al.ll2Qlll ) toge ther with the to- 
tal stellar mass, M* = (3±l)xlO lo M (|Swinbank et al.l 
12010ft . indicates baryon-dominated dynamics within the 
extent of the gas disk, (M gas +M*)/M dyn =l.l±0.2, and 
a gas fraction of / gas ~ M gas /M dyn ~ Af gas /(M* + 
M gas ) =0.6±0.1. This suggests a total baryonic mass 
surface density (stars and gas) within the disk of (3.0 ± 
0.5) x 10 9 M o kpc- 2 . 

The ratio of rotation/dispersion within a disk provides 
an important constraint on the amount of turbulence 
within the ISM. Nebular emission line studies have shown 
that high velocity dispersions are a common feature of 
high redshift galaxies (~4-10x higher than compara- 
bly luminous local galaxies) , possibly relat ed to the high 
gas fractions and star-formation rat e s fe.g.lLehnert et al.1 
12001 IForster Schreiber et al.|[200l iGenzel et al.ll201lD . 
However, measuring the typical local velocity dispersion 
is not trivial, especially with low (^ kpc) spatial resolu- 
tion, if the velocity gradient within the beam is much 
larger than the local velocity dispersion. In this case, 
beam smearing can produce a large bias in the recovered 
value (Davies et al. 2011). 

Since the spatial resolution obtained by our observa- 
tions is much higher than typically acheived in high- 
redshift galaxy studies and the orientation of the glaaxy 



is such that the maximum amplification is aligned 
roughly parallel to the velocity gradient, beam smearing 
has less of an effect on our intrinsic, although a small cor- 
rection is still needed. The maximum velocity gradient in 
SMM J2135-0102 is dV/dr = 0.15 km s" 1 per 1 , thus over 
a beam size of ~200pc there is a maximum contribution 
of 25kms _1 to the velocity dispersion. To account for 
this, we therefore correct the local velocity dispersion in 
each pixel and determine an average line of sight velocity 
dispersion of 91±12kms _1 across the galaxy image (we 
note that if instead we simply measure the uncorrected 
average line-of-sight velocity dispersion we derive a = 
97 ± lOkrns -1 ). Thus, we derive an average inclination- 
corrected v/a — 3.5±0.2. This is larger than that found 
in many of the star-form ing galaxies studied through 
their nebular emission (e.g. IForster Schreiber et aLll2009t 
lLaw et al.ll2009fl, and similar to tha t found in some local 
ULIRGs (fpownes fc SolonionlllQQSh . 

As an independent check of the kinematics of the 
gas reservoir in SMMJ2135 we use our CO(1-0) map 
from the EVLA. The resolution of this map is slightly 
poorer, but still comparable to that of the PdBI 
CO(6-5) observations, but importantly the CO(1-0) 
line emission encompasses much cooler and less dense 

gas (n^. 1 ^ ^400 cm -3 Ei/ke=5.5K) associated with non- 
star-forming ISM in addition to the warmer and denser 
gas, typically found in star-forming regions which is 

traced by the CO (6-5) emission (n^ ~6 x 10 4 cm~ 3 , 
E 6 /kB^116K). The velocity field and rotation curves 
from CO(1-0) in Fig. [2] show good agreement with 
those from CO (6-5) although it is clear that the 
CO(1-0) emission extends at least 1" (~ lkpc) fur- 
ther from the critical curve in both images of the 
lensed galaxy (even accounting for the slightly larger 
beam). This bias: much broader and more spa- 
tially extended low-J line emission compared to high- 
J lines has recently been reported for a number of 
high-redshift sub-mil limetre galaxies ([Ivison "eT!II[27)Tll : 
iRiechers et alJl2011allbTl in addition to low redshift star- 
bursts (e.g. iPapadopoulos fc Seaauistl 119981: iMao et al.l 
[20001: IWeifi et al.ll200lL 120051: IWalter et al.l 12002ft . A de- 
tailed comparison between low- and high- J CO emission 
and the radio continuum emission in SMM J2135 on sub 
100 pc scales from A-configuration EVLA observations 
will be discussed in an upcoming paper (R. Ivison et al., 
in prep). 

3.3. The high-pressure molecular ISM in SMM 32135 

The pressure in the ISM - overwhelmingly due to the 
non-thermal, macroscopic motions in the highly super- 
sonic, turbulent molecular gas - is an important pa- 
rameter for understanding the formation of stars. It 
effects the normalisation of the M(H2)-L C q( 1 _ q) rela- 
tion (iBrvant fc Scov illc 1996), Larson's relations (|Chiezd 
Il987t lElmegreenl Il989af) "and the star-formation effi - 
ciency A(SFR)/Am gas (ISchave fc Dalla Vecchial 120081) . 
It also plays a fundamental role in the Hi — > H2 
phase transition and thus in the formation of molecular 
cloud s in spiral disks (e.g. lElmegreenl[l993tlHonma et al.l 
19951) . and their fragmentation i nto dense cores (e.g. 
Chieze fc Pineau Pes Foretsl 11987ft . Finally, high pres- 
sure in a molecular inter-cloud medium may be responsi- 
ble for the shock compression of individual GMCs as they 
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fall into the central regions of ULIRGs and their s ubse- 
quent burst of star formation (|Solomon et al.lll997D . 

Before computing the average gas pressure in the ISM 
in SMMJ2135 it is important to note that the exis- 
tence of random macroscopic gas motions does not au- 
tomatically imply that the corresponding pressure term 
l/3p gas (cr 2 ) is physical and thus dynamically important. 
This requires a mechanism for momentum transfer be- 
tween the components to transmit the pressure, and one 
capable of doing so over a short "collision" mean-free 
path length with respect to the size of the gas distri- 
bution. Direct collisions between kinematically distinct 
components may not easily satisfy these two criteria in 
highly supersonic molecular clouds. On the other hand, 
magnetic field lines "threaded" through the various gas 
clumps may provide a suitable coupling between molec- 
ular clouds ([Elmegreenl U989bT ). However, it must be 
noted that the molecular gas in compact starbursts (e.g. 
ULIRGs - perhaps the closest analogues to SMMJ2135) 
may be very different from that found in the disk of the 
Milky Way, with a much larger - perhaps dominant - 
gas-mass fraction at d ensities of rt(H 2 ) > 10 5 cm -3 (e.g. 
IGao fe Solomonl 120041 ) . At such high densities and in 
the high-extinction environment of ULIRGs and SMGs, 
the average ionisation fraction of the gas may be very 
low, allowing ambipolar diffusion to very effectively re- 
move magnetic field lines from the gas, thereby reducing 
the coupling between components. Nevertheless, for the 
purposes of this study we assume that the computed gas 
pressure is physical and dynamically important. 

The dynamics and surface density of gas and stars can 
be used to estimate the ISM pressure, which can then be 
compared to that seen in gas-rich galaxies locally. For 
a rotating disk of gas and stars the external hydrostatic 
pressure at mid-plane is given by: 



Ptot ~ 



Egas + 









(1) 



where £ gas and £* are the surface density of the gas and 
stars, and tj gas and cr* are the vertical velocity dispersion 
of the gas and the stars, respectively (with no assump- 
tion that the gas scale leng th is smaller than the stars; c.f. 
iBlitz fe Rosolowskvl2006f ). F or the Milky Way this yields 
Ptot/ks ~ 1.4 x 10 4 cm" 3 K (|ElmegreenlH"989bl) . To esti- 
mate the stellar surface de nsity in SMM J2135, we use the 
stellar mass estimate from lSwinbank et al.l (|2010f) . M* = 
(3±l)xl0 10 M Q to derive £* = (2.4±0.5) x 10 3 M pc" 2 
(and assume that the stars and gas have the same spa- 
tial distribution, which is reasonable giv en the HST and 
Svitz er IRAC morphologies; Fig. Q] and Swinban k et al.1 
12010( 1 , whilst the gas surface density we derived above as 
E gas = 3.2 x 10 3 M Q pc~ 2 . 



The average velocity dispersion of the gas in the disk is 
Ogas = 91±10kms _1 , and we have to assume that the ve- 
locity dispersion of the stars is comparable, (although we 
allow this to vary between 50-150 km s _1 in the following 
results). The resulting mid-plane hydrostatic pressure is 
Ptot/k B ~ (2±l)xl0 7 Kcm- 3 . This pressure is - lOOOx 
higher than the typical ISM pressure in the Milky Way 
(10 4 Kcm -3 ), and approximately 2-3 x higher than seen 
in other local, gas-rich enviro nments, such as the Galac- 
tic Center or the Antennae ([Ro solowskv & Blitz] 120051 : 
iKeto fe MversllT986l : iWilson et alJl2003D . Only the com- 
pact disks in local ULIRGs have comparable mid-plane 
pressures due to th eir high gas-mass surface densities 
(> 5 x 10 3 M Q pc" 2 , IDownes fe Solomonl Il998h . 

We can compare our measured ISM pressures to 
estimates fo r five high-redshift disks galaxies from 
iGenzel et al.l ()2008l ) with well-constrained kinematics. 
Although gas masses do not exist for these galaxies 
directly, they have been inferred indirectly using the 
Ha luminosity density and the Schmidt-Kennicutt re- 
lation. This provides an estimate of the likely gas sur- 
face density to within a factor ~ 3x depen ding on the 
choice of parameters (e.g.fGenzel et al.ll2010l but see also 
Papadop oulos fe Pelupessvll2010l ). 

The (inferred) median gas mass for the IGenzel et al.1 
(2008) sample is M gas = 0.2xlO n M Q , with stellar 
masses of M* ~ (0.3-6) x 10 11 M Q in disks of scale 
length 4 -7 kpc and vertical scale heights of h* = 0.8- 
1.6 kpc ([Shapiro et al.l 120101) . Together, these proper- 
ties suggest ISM pressures in the range Ptot/fe ~ (0.3- 
10) x 10 7 Kcm -3 , suggesting that high gas pressures may 
be a ubiquitous feature of high-redshift, gas-rich systems. 

3.4. Disk stability and fragmentation scale 

Having establish that the ISM pressure is high, we next 
examine the susceptibility of the gas disk to collapse. In a 
rotating disk, perturbations smaller than a critical wave- 
length, L max , are stabilised against gravity by the veloc- 
ity dispersion; those larger than some critical wavelength, 
L m i n , are stabilised by centrifugal force. If the disper- 
sion and rotation velocity are too low, L m { n > L max and 
perturbations of intermediate wavelength grow exponen- 
tially. The Toomre parameter, Q = L max / L m i n , char- 
acterises the stability against local axisymmetric pertur- 
bations of a disk supported by differential rotation and 
random motion ([Toomrd fl964) . Galaxies with Q < 1 
are unstable on scales between L max and L min and will 
fragment into giant, dense clumps. This could trigger 
star formation in clouds of much higher mass and size 
than GMCs seen in local spiral galaxies with Q > 1. 
Dynamical friction, viscosity and tidal interactions may 
subsequently cause these clumps to migrate toward the 
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radius (kpc) 

Fig. 4. — Top: Two-dimensional kinematics of the molecular 
gas in SMMJ2135 from the CO(6-5) emission line after correct- 
ing for lensing amplification. The velocity field shows a veloc- 
ity gradient of ~ 500kms — 1 across ~3kpc in projection, indi- 
cating a dynamical mass of (6.0 ±0.5) X 10 10 Mq (corrected for 
inclination, i =60 ±8°). The contours denote the best-fit disk 
model and the crosses show the locations of the star-forming 
regions seen in the re st-frame 260-/mi continuum imaging from 
ISwinbank et all {2010). The source-plane beam is shown in the 
top right-hand corner and is approximately 150 X 850 pc on av- 
erage across the galaxy image, but reaches ~ 90 pc in the most 
highly amplified region. Inset: Velocity residual after subtracting 
the best-fit disk model. Bottom: Source-plane, one-dimensional 
velocity- and velocity-dispersion- profile from the CO (6— 5) emis- 
sion line extracted from the major kinematic axis using a synthetic 
slit of width 1 kpc. The dashed curve denotes the rotation curve 
of the best-fit two-dimensional velocity field. The solid horizon- 
tal bars denote the resolution of the CO (6— 5) observations as a 
function of position within the galaxy. The lower panel shows the 
one-dimensional velocity dispersion across the same major kine- 
matic axis of the galaxy (small points) whilst the the large solid 
points show the velocity dispersion of the individual star-forming 
regions as identified from the SMA morphology. 



center of the galaxy potential, forming a bulge which 
stabilises the syste m against further fragmentation (e.g. 
IGenzel et al.lf2008h . 

The Toomre parameter is calculated from 



where rj r is the one-dimensional random velocity disper- 
sion, k is the epicyclic freq uency and E is the surface 
mass density (jToomrd fl96l . The value of k is uncer- 
tain as it depends on the unknown mass distribution; 
our observations are consistent with a range V^r - 
corresponding to constant V c and V c oc R, respectively. 
Disk thickness and gas abundance also affect the value 
of Q, and we have assumed that the measured velocity 
dispersion is equal to a T . Adopting k = V3 V "r* ap- 
propriate for a uniform disk, and using the dynamical 
mass to estimate the surface mass density, E, we find an 
inclination-corrected Q = 0.50 ± 0.15 (the error is con- 
tributed equally from the uncertainty in the dynamical 
mass and the assumed k). The gas disk in SMMJ2135 
thus appears to be dynamically unstable (Q < 1), per- 
haps more gravitationally unstable than the compact, 
gas-ric h disks residing in ULIRG s where Q is closer to 
unity (jDownes fc Solomon|[l998l) . 

Empirically, this instability will cause large, dense gas 
condensations to form in the molecular gas distribu- 
tion. This fragmentation should occur on the scale of the 
Jeans length for dispersion support. In a uniform disk, 
the largest scale for which velocity dispersion stabilises 
against gravitational collapse is: 



(3) 



8GE 



ttGE 



(2) 



which can be estimated from the dispersion and dynami- 
cal mass density. As with Q, the unknown mass distribu- 
tion, disk height, gas content and directional dependence 
of a result in an uncertainty of around 30 per cent. The 
resulting instability scale is Lj = 400 ± 150 pc and the 
corresponding unstable mass scale is cr 4 /G 2 E ~ 1O 9 M , 
equivalent to the entire molecular gas reservoir of a spiral 
disk such as the Milky Way. 

Our previous high-resolution, rest-frame far-infrared 
imaging has shown that the galaxy contains several 
bright star-forming regions, with star-formation rates 
of 30-100 M Q yr _1 each (equivalcntly, luminosities of 
2-5 x 10 11 L(T ) ), spre ad across a radius of ~ 2 kpc 
(|Swinbank et al.ll201Q[) . Although the velocity integrated 
CO (6-5) emission line maps appear smooth, Fig.[3]shows 
that when the cube is viewed in the three dimensions 
the CO(6-5) and CO(1-0) have complex morphologies. 
Hence, to examine the structures in the gas distribu- 
tion and so examine how these might relate to the 
rest-frame 260/im continuum morphology, we extract 
position- velocity plots and the one-dimensional intensity 
profiles of the CO(6— 5) and CO(1-0) emission from these 
PV diagrams and also show these in Fig. G3 These show 
that the gas morphology is indeed structured, apparently 
c omprising (at lea s t) five clumps across the ^5 kpc disk. 

iDanielson et al.l (|2011f) demonstrated that to model 
the integrated multi-transition CO kinematics of 
SMM J2135 required multiple velocity components which 
they denoted X, Y, Zi and Z 2 . Our high-resolution CO 
dynamics now allow us to spatially locate these struc- 
tures and link them to the star-forming regions seen in 
the rest-frame 260-/Ltm emis sion. Fig. [2]fc[3]show th at the 
lowest velocity material in Daniclso n et al.l ()201lD ."Zi" 
and "Z 2 " (at Au=28±9 and -167±9kms _1 with re- 
spect to the systemic respectively) appear to relate to 
two clumps separated by < 0.2" in projection on the 
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Fig. 5. — The CO(6-5) and CO(l-O) spectra extracted at the positions of the star-forming regions identified by the peaks in the rest-frame 
260-/xm emission. The datacube has first been corrected for the large-scale velocity field in order to remove the broadenin g of the lines due 
to large scale motions. Components Z1Z2, Yi, Y2 and "X" correspond to kinematic components "Z" and "Y" and "X" in Daniclson et al. 
(2011) (and SMA clumps "Z, Y & X" in Fig. |2j. The clumps have velocity dispersions of 45-85 kms' 1 which is a factor ~ 5x higher than 
expected for local, quiescent GMCs given their size and luminosity. We also overplot the equivalent CO(1-0) line widths, which broadly 
trace the same kinematics as the CO (6— 5) emission. 



sky, are associated with the rest- frame 260- ^.m emis- 
sion cl osest to the criti c al cur ve. Kinematic component 
"Y" in lDanielson et al.l (|2011[ ) (Aw=165±13kms- 1 ) ap- 
pears to comprise two or three gas rich clumps between 
v ~ 1 10-280 km s -1 . This gas-rich component is centered 
within ~ 0.2" of the brightest emission at rest- frame 
260 /xm. Finally, the highest velocity material, "X" in 
IDanielson et al.l (| 2 1 If) (Av=396±9kms~ 1 ) corresponds 
to the highest velocity, faintest emission in the rest-frame 
260-/xm and CO(1-0) and CO(6-5) maps. 

Given the correspondence between peaks in the molec- 
ular gas and the rest-frame 260-/xm emission, we will use 
the rest-frame 260-/xm map as a guide to isolate the star- 
forming regions, extracting the properties of the molecu- 
lar gas at these locations. Clearly this biases our results 
towards the properties of the brightest star-forming re- 
gions, potentially missing some of the gas emission lies 
outside these actively star-forming regions, however, it 
has the advantage that we do not need to identify peaks 
directly from the CO emission - which itself is difficult 
to accomplish, without introducing biases due to surface 
brightness effects and signal-to-noise. We therefore re- 
strict the following analysis to the locations of the bright- 
est regions in the rest-frame 260- /an emission. 

Before extracting the size, luminosity and velocity dis- 
persion of the molecular gas we must first remove the 
large-scale velocity structure of the galaxy (and hence 
any rotational contribution to the line widths of the 
clumps) , as well as estimating the contribution of the ve- 
locity gradient within the beam. We therefore subtract 
the best-fit kinematic model for the galaxy velocity field 
from the cubes and then extract CO(6-5) and CO(1-0) 
spectra and measure the line luminosity, velocity disper- 
sion and central velocity (Fig. [SJ. As in § I3.2l we use the 
model velocity field to estimate the local velocity gradi- 
ent (dV/ dR) within the beam and remove this in quadra- 
ture from the velocity dispersion. The largest corrections 
are in the central regions of the galaxy where the local ve- 
locity gradient is large (25kms _1 over 200 pc). We note 
that in all four regions from where we extract spectra, 
the CO has intrinsic velocity dispersions of 45-82 km s _1 
from CO(6-5), and that the CO(1-0) line widths are a 
factor 1.2 ±0.1 times higher (even after beam-matching 
the cubes). However, this is as expected if the CO (6- 
5) traces denser gas regions, more deeply embedded in- 
side molecular clouds, and thus subjected to higher av- 
erage pressures. To estimate the approximate sizes for 



the molecular gas associated with the star-forming re- 
gions, we collapse the CO (6-5) source-plane datacube 
over the ±2xFWHM of the CO line at its systemic ve- 
locity and fit an elliptical Gaussian profile to the result- 
ing image - we use the CO (6-5) here since the resolution 
is slightly higher than the CO(1-0). The resulting sizes 
in the source plane are rco= 1 10-300 pc, comparable to 
the Jeans length computed in §3.3 for the gaseous ISM in 
this system, which would be the expected from the initial 
fragmentation scale of a gravitationally-unstable gaseous 
disk. Since the source-plane reconstruction is sensitive to 
the lensing model, we also reconstruct and measure the 
sizes of each of the star-forming regions from the family 
of lens models which adequately reproduce the lensing 
configuration, which, on average introduces a ~50% er- 
ror on this size measurement. Since measuring the sizes 
for gas emission is difficult at these surface brightness 
levels, for all size measurements, we conservatively allow 
a factor 2x error, and report these in Table 2. Finally, 
we estimate the molecular gas mass in the vicinity of the 
star-forming regions using the CO(1-0) line luminosity, 
obtaining M = (4-15)xl0 8 M Q for a=2. Together, the 
gas in the vicinity of the brightest star-forming regions 
therefore comprises approximately 10-20% of the total 
gas in the system. 

Finally, before investigating the properties of the star- 
forming regions in detail, it is useful to compare their 
characteristics with those in comparably luminous galax- 
ies where comparab le resolution observa t ions h ave been 
made. In particular, Downes fc Solomonl (|1998t ) (see also 



Saka moto et al.|[2008 ) show that the molecular gas within 
the central regions of local ULIRGs are characterised by 
Q < 1, typically with two to three compact, dense re- 
gions with radii of ~ 70-100 pc and masses ~ 10 9 M© 
(and hence densities, ~ 2x 10 4 cm" 3 ) and typical infrared 
luminosities of (0.3-1) x 10 12 L Q . These properties are 
similar to those for the clumps in SMMJ2135. However, 
the key difference is that the bright star-forming regions 
in local ULIRGs tend to be located within the central 
50-200-pc radius of the disk, whereas the star-forming re- 
gions within SMMJ2135 are distributed across ~ 2kpc. 
This shows that while the star formation is occuring in 
similarly compact clumps, these clumps are more widely 
distributed across a much larger gas disk in this high- 
redshift galaxy. 
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3.5. Cloud scaling relations in a high-pressure 
environment 

Next, we compare the size, luminosity and velocity dis- 
persion of the star-forming regions to local GMCs and 
starburst complexes. For the line luminosities and ve- 
locity dispersions, we use the average from the CO(l- 
0) and CO (6-5) measurements, but ensure that the er- 
rors include the range of solutions from each of the two 
transitions separately. In Fig. [6] we compare the size 
and velocity dispersion of the molecular gas associated 
with the star-forming regions. This figure shows that the 
star-forming regions do not lie on the local relation fo r 
GMCs in quiescent environments (|Bolatto et al.l 12008). 
but are instead systematically offset such that at a fixed 
size, the velocity dispersion is ~ 10 x higher than typi- 
cal GMCs. Moreover, the offsets of the clumps seen in 
Fig. [JJ] from local scaling relations is comparable to that 
found for GMCs in the Galactic Center and other gas- 
rich environments. These offsets have been interpreted 
as a consequence of the high external ISM pressure on 
the c l oud s urfaces (jKeto et alJl2005b iBlitz fe Rosolowskyl 
l27Ml200l . 

We can compare the offsets in the size-velocity rela- 
tion with theoretical expectations for the line width-size 
relation for GMCs in high-pressure ISM environments. 
Two important scaling relations deduced for molecular 
cloud s in the Gala xy are the so-called Larson's scaling 
laws (|Larsonl [1981). namely: a) the CO velocity line- 
width-size relation a = AR a , and b) the average gas 
density-size relation (n(H2)) = BRr k . The first is the 
most robust of the two, and has been found to be remark- 
ably constant in the Galaxy, both in terms of its n ormal- 
isation and power-law index (jHever fc Brunt 2004) . This 
reflects the near-universality of turbulence in molecular 
clouds, while suggesting that the energy injection neces- 
sary to drive the turbulence occurs mostly at low spatial 
frequencies (i.e. large scales). The turbulent molecular 
clouds found to obey the Larson re lations are almost al- 
ways near-virialised (|Larsoniri981l ): in such a case, the 
density-size power law can be readily derived from the 
line-widt h-size relat i on and the virial theorem. In deed 
following lElmegreenl (|1989bD (see also lChieze]|1987ft . ap- 
plication of the virial theorem on molecular gas regions 
of radius, P, and boundary pressure, P ex t yields 



-Pcxt/fcB 

10 4 KcnT 



1/4 



1/2 



and 



(n(H 2 )) = n 



10 4 K cm" 



1/2 



(4) 



(5) 



where a Q = 1.2 km s 1 and n ~ 10 3 cm 3 , as o btained 
from studies of molecular clouds in the Galaxy (|Larsonl 
H98l IWolfire et al.l [20031: IHever fc Brunt! l200l . Theo- 
retical studies predict similar values t o within a factor 
of two (jChiezd .1987: El megreenlll989bj). Equatio n 151 can 
be recast as a mass-radius relation (jElmegreenl 11989b) 
and derived from the velocity-line width relation and the 
virial theorem applied for gaseous regions with boundary 
pressure. For the chosen normalisation value of the for- 
mer, n , this yields: 



M(H 2 ) = 290 



10 4 C 



; K 



1/2 / R x 2 

-) M . 

pcy 



(6) 



The large mid-plane hydrostatic pressure, Ptot/^B ~ 
(1 - 3) x 10 7 Kcm- 3 in SMMJ2135 will correspond to 
an external cloud boundary pressure P cx t ~ Ptot/(l + 
o-o + (So) = (0.6 — 1.8) x 10 7 cm~ 3 K for relative cosmic 
ray and ma gnetic pressure con tributions of a = 0.4 and 
f3 Q = 0.25 (Flmcgrccn 19891], and references therein). 
From equations H& [5] we then expect a (SMMJ2135) ~ 
7 x cr (Galactic) and n (SMMJ2135) ~ 50 xn (Galactic). 
The former is in very good agreement with the offset for 
the star-forming regions (a factor 7-10 offset in a at a 
fixed size), whilst the latter implies that the molecular 
gas is expected to be ~ 50 x denser on all scales when 
compared to the gaseous ISM in the Milky Way. Thus 
for a typical GMC size of 2P ~ 20 pc in SMMJ2135, 
the average densities would be ~ 5000 cm -3 rather than 
— 100 cm -3 as in the Milky Way. 

In Fig. [75] we also plot the mass-radius relation for a 
range of turbulent ISM pressures and overlay observed 
values of M and R from Galactic GMCs, ranging from 
those in the disk to those in the high-pressure environ- 
ment of the Galactic Center. The star-forming molecu- 
lar gas regions of SMM J2135 are offset by two orders of 
magnitude with respect to quiescent GMCs, as expected 
from Equation [6j and the turbulent pressures deduced 
for its gaseous ISM. 

Whilst our analysis could readily attribute the off- 
sets of the line width-size relation to the independently- 
derived large ISM pressures, we caution that we may 
still be probing scales at or just beyond where this re- 
lation is valid. Indeed this scaling law is expected to 
have a cut-off at the scales corresponding to the largest 
possible virialised molecular cloud structures, or simply 
to the lowest spatial frequency of turbulent driving (the 
disk scale- hei ght being an obvious choice). The origi- 
nal study by iLarsonl (|1981l) extends up to a cloud size 
of 100 pc, while the mos t recent and systematic study 
by IHever fc Brunt! (|2004l ) only goes up to 50 pc (a typ- 
ical GMC size). On larger scales, systematic (stream- 
ing) motions between otherwise virialised gas structures 
could broaden the apparent "cloud" line widths, our re- 
moval of regular disk-like motions notwithstanding. Such 
systematically-broadened line widths could then create 
apparent offsets from the local scaling relations, simi- 
lar to those seen in Fig. HI that are unrelated to in- 
creased ISM pressures. Only ALMA will be capable 
of verifying the true pressure-induced offsets in distant 
star-forming galaxies, by probing scales <100pc with 
sufficient sensitivity and spatial resolution to establish 
both the slope and the normalisation of the ISM scaling 
laws at high redshifts. This will be possible for strongly- 
lensed starburst systems such SMMJ2135, allowing an 
unprecedented insight into key quantities characterising 
the turbulent molecular gas in star-forming systems at 
high redshifts. 

3.6. Star-formation efficiencies at high ISM pressures 

High average volume densities for molecular gas have 
been revealed by HCN, CS and HCO + observations 
of local ULIRGs and interpreted as being responsible 
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Table 2: CO properties of the clumps 





Zi 


z 2 


Yi 


Y 2 


X 


Galaxy Integrated 


S , 870Aim (mJy) 


8.2±2.1 


8.2±2.1 


6.4±2.1 


19.1 ±2.1 


7.5 ±2.1 


106.0 ±3.5 


Amplification (fi) 


63 ±14 


63 ±14 


53 ±5 


40±4 


20 ±2 


37.5 ±4.5 


SFR (M yr" 1 ) 


31 ±10 


31 ±10 


30 ±10 


97 ±14 


70 ±14 


400 


Vel (kms" 1 ) (z sys = 2.32591) 


-185 ±19 


-63 ±9 


116 ±12 


280 ± 20 


380±20 




Sco(l-O) (Jykms" 1 ) 


0.06 ±0.01 


0.05 ±0.01 


0.08 ±0.01 


0.11±0.01 


0.14 ±0.01 


2.12 ±0.02 


L CO(l-0) (xl0 8 Kkms- 1 pc 2 ) 


2.0±0.6 


1.7±0.6 


3.9±0.4 


7.1 ±0.8 


1.8 ±0.2 


173 ±9 


^00(6-5) (Jykms" 1 ) 


0.57 ±0.05 


1.40 ±0.05 


1.53 ±0.05 


1.80 ±0.05 


1.02 ±0.05 


21.5±1.1 


L CO(6-5) (xl0 8 Kkms- 1 pc 2 ) 


0.54 ±0.05 


1.30 ±0.06 


1.94 ±0.07 


3.1 ±0.09 


3.41 ±0.17 


48 ±2 


°"CO(6-5) (kms" 1 ) 


45 ±15 


82 ±7 


63 ±12 


74 ±12 


56±15 


210 ±30 


CT CO(i-0) (kms" 1 ) 


47 ±20 


110 ±20 


57 ±14 


79 ±15 


78±15 


210 ±30 


r CO(6-5) (PC)* 


110 


116 


105 


160 


300 




rs70)im (pc) 


90 ±20 


90 ±20 


98 ±30 


180 ±230 


[352,452] ± 50 




M gasCO (i-0) (xlO 8 M ) 


4.0±1.5 


3.3±1.5 


7.8 ±0.8 


14.2 ±1.6 


3.6 ±0.4 




M gasCO (6-5) (xlO 8 M ) 


3.8±0.4 


9.2 ±0.8 


13.8 ±1.0 


22.8 ±1.0 


24.3 ±1.6 





TABLE 2 

Notes: the line luminosities, L', are corrected for lensing ampl ification. The gas mass (M ga a) is calculated from 
Mgaa = Q ^cO(i-0) WITH Q = 2 (and Rgi = 0.28 ± 0.03 for the CO(6-5): IDanielson et al"]|2011I) , For ROWS marked with *, we 

estimate that the uncertainty on the measurement is a factor ~2. 



for the very high star-formation efficienc ies observed 
for t heir global molecular gas reservoirs reve et aLl 
2009). The high luminosities of the star- forming regions 
in SMMJ2135 per gas mass suggest high efficiencies, 
with EFiR/Sg as = (100-450) L Q /M Q , near the radiation- 
pressure limit recently advoc ated for extreme starb ursts 
(|Andrews fc Th ompson 2011]; see also lScovillell2004l for a 
first derivation). 

However, with estimates for the size, velocity disper- 
sion, mass and star-formation rates of the clumps in 
SMMJ2135, we can estimate the star-formation efficien- 
cies per dynamical time directly. Using the radii of 
B ~ 100-200pc (as measured from the CO(6-5) line 
emission and the rest- frame 260-/iin emission), for ve- 
locity dispersions of 40-60 kms -1 , the crossing time of 
the clump region is r ~ 4-5 Myr. The inferred star- 
formation rates of each region is ~ 30-90 M Q yr -1 , thus 
in one dynamical time, l-5xlO 8 M0 of gas will be con- 
verted into stars. Assuming that all of the molecular gas 
in a GMC is converted into stars (i.e. no gas loss), it 
takes 2-10 dynamical times to convert 10 9 Mq into stars. 
Conversely, the star-formation efficiency per dynamical 
time must be e c i ~ 10-50 per cent, i.e. approaching the 
star-formation efficiency of dense HCN-bright cores in 
the Milky Way. This is indeed a natural outcome of the 
much higher densities expected at every spatial scale in 
the hig h-pressure ISM within SMM J2135 (as discussed 
in £I3.5|) . Clearly these calculations are simplistic since 
they assume that there is no gas loss and that the star- 
formation rate remains constant over the gas depletion 
lifetime. Nevertheless, they demonstrate that the star- 
formation efficiency in the clumps can be much higher 
than those in molecular clouds in quiescent GMCs. It 
is worth noting that the galaxy-averaged star-formation 
efficiency is much lower, e gat ~ e ct x M c i umps / M to tal ~ 2- 
10% (i.e. closer to t he average value seen in the Schmidt- 
Kennicutt relation; lKennicuttlll998l) . 

We can compare these values to a theoretically-derived 
maximum efficiency expected for high-pressure molecu- 
lar ga s in star-forming regions. lElmegreen fc Efremovl 
(1997) suggest that if both the gas density and the rate 



of star formation is high enough, most of the gas in 
a cloud may be consumed before th e mole cular cloud 
is dispersed. Elmcgr e'en fe Efremovl (|1997t ) also show 
that high-pressure environments result in a lower dis- 
persal rate, which improves the efficiency (resulting from 
the higher physical density and gravitational binding en- 
ergy). In this model, the star- formation efficiency de- 
pends on the relative rate of star formation (which is 
proportional to some power of the gas density, n ~ 1.5; 
lKennicuttlll998l ) and the rate of dispersal (which is pro- 
portional to the stellar luminosity and inversely propor- 
tional to the gravitational binding energy). By scaling 
the relative rates of star formation and efficiency from 
those in the Milky Way, lElmegreen fc Efremovl (|1997f ) 
estimate the star-formation efficiency as a function of 
pressure and cloud mass, and in Fig. [7] we show how the 
star-formation efficiency scales with ISM pressure and 
mass. For typical GMC mass scales of £ 10 5 M Q and 
external ISM pressures of 10 3 x that of the Milky Way, 
the star-formation efficiency can be as high as ~ 80 per 
cent. 

4. IMPLICATIONS FOR STAR-FORMATION IN HIGHLY 
TURBULENT GAS DISKS 

The large line width observed for the high-pressure 
molecular gas of extreme starbursts like SMM J2135 will 
alter the so-called sonic or turbulent pressure length 
A s . This marks the scale at which supersonic turbu- 
lence becomes subsonic and near-thermal, and below 
which micro-turbulent velo c ity fields dominate the ga s 
motion (jWolfire et all 120031 : iKrumholz fe McKed l2005h . 
In other words A s marks the smallest scale where the 
line width-size relation can hold in molecular clouds, be- 
low which there are near-constant velocity dispersion s 
that no longer correlate with sizes (PI umc 
Using a line width-size relation re-cast in terms of the 

1 /2 

isothermal sound speed c s , namely <r(R) — c s (R/X s ) , 
it is clear that a ~ 10 x higher normalisation factor 
implies a difference in A s between SMMJ2135 and the 
Milky Way of Af MM /A^ w ~ 0.01 (assuming the same 
sound speed). This means that in systems such as 
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Fig. 6. — The molecular cloud scaling relations for the star-forming regions within SMMJ2135 compared to those in galaxies in the 
local Universe. Left: The velocity dispersion versus cloud radius. The extragalactic GMCs in quiescent environments define the local line 
width-size relation (Larson 198f ), with the normalisation consistent with the average mid-plane pressure in the mid-plane of the Milky Way 
(Elmegreen 1989). Those GMCs in gas-rich, high-turbulent-pressure environments (such as the Galactic Center) tend to be systematically 
offset from this relation. The line width-size data for star-forming regions in SMM J2135 are compatible with a line width-size relation 
with a ~ 5x higher normalisation than found in the Milky Way. Right: The M(H2) — R 2 relation (equivalent to the second cloud scaling 
relation expressed in Eq. [5}. The molecular star-forming regions in SMMJ2135 lie on the line corresponding to the very high turbulent 
pressures estimated for them (Equation [TJ. 



SMM J2 135, much smaller gaseous structures can re- 
main supported against gravitationally-induced fragmen- 
tation/collapse by strong turbulence is dominant on ~ 
100 x smaller scales than in the kinematically quiescent 
ISM of the Galaxy. 

In the context of turbulence-regulated star formation, 
only regions with overdensities x = n/(n) can "escape" 
the grip of supersonic turbulence and irreversibly pro- 
ceed towards gravitational collapse ([Krumholz fc McKeel 
2( )()•">;•. For the Milky Way, the critical value is 
Xcrit = 1.07M 2 (where M is the one-dimensional 
Mach number). For SMMJ2135, our results suggest, 
x cr it (SMM J2135)~100x cr it(Galaxy). I n absolute den- 
sities , for (n)(SMM J2135)-10 4 1 cm" 3 (jDanielson et al.l 
120111 ). the critical density will be n crit ,gr~l-4 x 10 8 cm~ 3 . 
This is an exceptionaly high density, and for perspective 
the corresponding critical density in the Milky Way is 
~10 4 -10 5 cm -3 . As such, the normalized star-formation- 
rate ( SFR) per free fall time wi ll be increased substan- 
tially (|Padoan fc N ordlund 20111 ). Far from being the ex- 
ception, star formation in such highly turbulent ISM will 
not be atypical in the distant Universe where starbursts, 
such as SMM J2135, represent an important star-forming 
population. Indeed, half of the stars seen today m ay have 
formed in such environments (jBlain etaII[T999l) . 

4.1. The thermal state for dense highly turbulent gas 

The higher normalization of the linewidth-size relation 
in the ISM of SMM J2135 translates to a much smaller 
sonic length controling the onset of near-thermal motions 
and eventual gravitational collapse. However this as- 
sumes isothermal gas, and thus an invariant sound speed 
c s . While th is has been widely used in numerical simu- 
lations (e.g. IPadoan fc Nordlund! I2002D . and analytical 
models of turbulence-regulated star-formation in indi- 
vidual molecular clouds a nd gaseous galactic disks (e.g. 
iKrumholz fe M cKcc 2005|), it does not represent a real- 
istic ISM. Indeed, non-isothermal gas may play a crucial 



role in defining the thermodynamic state of the dense 
gas and hence on t he stellar IMF and its mass scale 
(i.e. the IMF "knee " ; ILi et al.|[200l Uappsen et al|[2005l: 
IBonnell et all[2006l1 s l 

The very strong turbulence present in the ISM of galax- 
ies such as SMMJ2135 (as well as local ULIRGs) will 
influence the thermodynamic state of their UV-shielded 
dense gas phase (where the initial conditions for star for- 
mation are set) in two important, yet opposing ways. 
First, it will volumetrically heat this gas phase (i.e. un- 
like the surface-like heating of far-UV, optical or infrared 
photons), and second it will enhance molecular line cool- 
ing in much deeper and denser regions of the ISM. The 
latter occurs simply because highly supersonic line widths 
can now be maintained in much smaller, and thus denser 
gas regions (cf. Eq. [5]). This will reduce the radiative 
trapping of cooling molecular line photons from such re- 
gions deep inside molecular clouds, as molecular line opti- 
cal depths, r cx (dV/dR)" 1 , can now remain low at much 
higher densities because of the much larger velocity gra- 
dients, dV/dR. These two effects can significantly alter 
the initial conditions for star formation in the highly tur- 
bulent ISM by changing the dense core mass spectrum, 
and/or the density, n c setting the equation of state inflec- 
tion point (the latter is the density where the polytropic 
index flips from <1 (for n < n c ) to > 1 ( for n > n c ), i ts 
value may determine the IMF mass scale (|Larso n 200j)). 

The effects of large velocity gradients on the degree 
of gas fragmentation and the resulting mass function of 



The polytropic Effective Equation of State Pthcrmal — 
is a simple way to incorporate the complexities of strongly evolv- 
ing gas heating and cooling processes in such numerical simulations 
and the values of its polytropic index 7 determine the degree of gas 
fragmentation and the resulting den se gas core mass spectrum (e.g. 
ILi et al.l l2003; Jappscn et al. 2005). The dependence of 7 on the 
properties of the ambient ISM (e.g. metallicity, velo city field, back- 
ground radiation) have been explored in some detail ISpaans fc Silkl 
I2000D . 
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Fig. 7. — Efficiency of star formation as a function of cloud mass 
and ISM pressure (Pq d enotes the Milky Way value). The lines 
show the calculation from Elmcerccn & Efrcmov ( 1997) for the ef- 
ficiency of star formation as a function of the cloud mass in an 
environment with pressure 10 3 X that of the Milky Way (the effi- 
ciency as a function of cloud mass for the Milky Way mid-plane 
pressure, log(P/PQ = 0) is also shown). For the ISM pressure mea- 
sured in SMM J2135 (log(P/P @ = 3) and typical GMC masses of 
^, 10 5 Mq , the star-formation efficiency can be as high as ~ 80-90 
per cent. 

dense cores have been explored, but only up to the rather 
low v alue of dV/dR ~ 3kms _1 pc~ 1 (Sp aans fc Silk! 
I2000D . In the case of SMM J2135, where strong turbu- 
lence can be maintained at n(H2) > 10 5 cm~ 3 , the aver- 
age velocity gradient expected from self-gravity-induced 
gas motions alone will be: 



dr 



n(H 2 ) 
10 3 cm" 



1/2 



6.5\/akms 1 pc 



and may approac h dV/dR 



ap proac h av/dH ~ 10kms _1 pc~ 1 (see also 

iGreve et al.l [2009). To date, such strongly turbulent 
media have not been si mulated, even for a purportedly 
starburst-like ISM (e.g. IKlessen et all 120071) . The dis- 
covery of the extreme kinematic properties of the ISM in 
SMM J2315 makes such simulations necessary, with pre- 
liminary results indicating significant implications for the 
resulting dense gas fragmentation and its star-formation 
efficiencies (Hocuk & Spaan, in prep). 

5. CONCLUSIONS 

We have mapped the distribution of molecular gas in 
the lensed, star-forming z — 2.3 galaxy, SMMJ2135, 
through the CO(6-5) and CO(1-0) emission lines, us- 
ing the PdBI and EVLA respectively. After correcting 
for lensing, our high-resolution data provide a physi- 
cal resolution of ~ 100-200 pc. The dynamics of both 
the CO(1-0) and CO(6-5) line-emitting material con- 
vincingly demonstrate for the first time in a galaxy at 
z> 1, that the molecular gas is located in a large rotat- 
ing disk, which suggests such structures may be ubiqui- 
tous in high-redshift galaxies. From the kinematics we 
derive an inclination-corrected rotation speed for the disk 



of v Iot = 320 ± 25kms 1 , v/a — 3.5 ± 0.2 and a dy- 
namical mass of (6.0 ± 0.5) x 10 10 M Q within a radius 
of 2.5 kpc (the spatial extent of the detected emission in 
CO(1-0)). The disk is massively unstable, with a toomre 
parameter, Q = 0.50 ±0.15, which is lower than the com- 
pact, gas-rich disk s residing in local ULIRG s where Q 
is closer to unity ([Downes & Solomon! [19981) . We also 
find the CO (1-0) emission is slightly more spatially ex- 
tended than than the higher- J CO lines, as has been seen 
for other high-redshift SM Gs and ULIRGs (jlvison et all 
l201lHRiechers et alJl20TlaMbh . Comb ining with measure- 
ments of the stellar mass within this system, we find 
that the dynamics of the disk are baryon dominated, 
(M gQS 4-M*)/M t ;j / „=l.l±0.2, with a molecular gas frac- 
tion of Mg as /(Mg as + M*) = 0.6 ± 0.1, which is compa- 
rable to that found in other starburst galaxies at t hese 
epochs (jTacconi et al.|[2008l [2oTot lEngel et al.ll2010D . 

We use the gas dynamics and surface mass density to 
explore the ISM properties, deriving a mid-plane hydro- 
static pressure of Ptot ~ (2 ± 1) x 10 K cm -3 , which cor- 
responds to an external boundary pressure on the GMCs 
of P oxt = (0.6-1.8) x 10 7 cm- 3 K, a factor -1000 times 
higher than the Milky Way. Using established expres- 
sions of the (linewidth)-scale and (average density)-scale 
scaling relations known in the local ISM we deduce that 
at these high pressures the expected velocity dispersions 
and densities of molecular gas regions should be 7 x and 
~ 50 x the Galactic values respectively at all scales where 
supersonic turbulence remains important. We use our 
high-resolution CO cubes to demonstrate that this is in- 
deed the case for the star- forming regions in SMM J2135, 
which are similar to the high-pressure ISM environments 
found in the Galactic Center. 

The systematically higher gas densities expected from 
the higher normalization of the (average density)- (scale) 
relation (or its equivalent (cloud mass)- (radius) relation) 
are expected to dramatically elevate the star-formation 
efficiencies in this system, providing a natural explana- 
i (7)tion for the high star-formation-rate densities within the 
compact regions seen in the high resolution rest-frame 
260/im continuum imaging of this galaxy. We isolate 
these star-forming regions using their rest-frame 260-^m 
continuum emission and extract the gas velocity disper- 
sions and luminosities in the corresponding regions of the 
disk. The molecular mass in the vicinity of each of the 
star-forming regions is M c \ = (5-15) x 1O 8 M . Thus, 
the brightest star-forming regions in the galaxy make up 
approximately — 10 per cent of the total baryonic mass 
in th e disk (see also lGenzel et alJl201fl : lElmegreen et al.l 
[2009h . In all four regions, the velocity dispersions of the 
molecular gas are 40-85 km s -1 . 

We then compare the measured star-formation efficien- 
cies of the clumps with those from theoretical expecta- 
tions of high pressure ISM. Given their sizes and veloc- 
ity dispersions, the crossing time of each of the star- 
forming regions is r — 4-5 Myr. For star-formation rates 
of — 30-90 M Q yr -1 , each star-forming region will con- 
vert the ~ 10 9 Mq of molecular gas in their vicinity into 
stars in just 2-10 dynamical times. Conversely the star- 
formation efficiency per dynamical time must be ~ 10-50 
per cent, a factor ;> 10 x higher than in the GMCs in the 
Milky Way. Nevertheless, this is compatible with the 
theoretical expectation for the star-formation efficiency 
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at hi gh-pressure and cloud mass ([Elmegreen fc Efremovl 
[1991 . 

Finally, recasting the linewidth-scale relation in terms 
of the sound speed in an isothermal medium, as well as 
in the context of turbulence-regulated star formation, we 
show that for the highly turbulent ISM of SMMJ2135 
the turbulence will remain dominant, preventing gravi- 
tational collapse for structures that are up to 100 times 
smaller (and denser) than the corresponding structures 
in the Milky Way, resulting in a high critical density of 
~10 8 cm -3 for the onset of star formation. In the context 
of a more realistic non-isothermal ISM, we argue that 
the highly supersonic gas motions maintained at much 
higher gas densities in the turbulent ISM may alter the 
thermal state and the effective equation of state of the 
gas phase from which stars form. Detailed numerical sim- 
ulations of such highly turbulent ISM are thus urgently 
needed to explore its mass fragmentation spectrum, the 
star-formation efficiencies, and the possible impact on 
the emergent stellar IMF. 

Overall, our observations have resolved the molecular 
gas emission in a high-redshift star-forming galaxy at 
z = 2.3 on 100-200 pc scales. The gas is located in a 
rotating disk which is highly unstable (Q=0.50±0.15), 
and the gas density suggests the ISM should be under 
a mid-plane hydrostatic ^10 3 x that of the Milky Way. 
We show that the gas in the vicinity of the star-forming 
regions is ~10x denser than that of typical GMCs in 
the Milky Way, and causes them to be significantly off- 



set from scaling relations which govern the structure of 
local GMCs. These offsets imply that supersonic tur- 
bulence will remain dominant on scales ~100x smaller 
than in the kinematically quiescent ISM of the Milky 
Way, while the molecular gas will be ~50x denser on 
all scales. Thus, the high star-formation densities seen 
in SMM J2135 are a consequence of the high gas surface 
densities and pressures, which result in a low dispersal 
rate for the clouds and star- formation efficiencies ~ 10 x 
greater than those in quiescent environments. In the era 
of ALMA under full science operations, such observa- 
tions should become common place, verifying the effects 
of pressure-induced offsets from local GMC scaling rela- 
tions in distant star-forming galaxies. 

ACKNOWLEDGMENTS 

We would like to thank the referee for a constructive 
report which significantly improved the content and clar- 
ity of this work. AMS, IRS, RJI and APT gratefully 
acknowledges support from STFC. We thank Richard 
Bower, Bruce Elmegreen, Natascha Forster-Schreiber, 
Reinhard Genzel, Phil Hopkins, Mark Krunholz, Norman 
Murray, and Linda Tacconi for advice and conversations, 
and Jim Geach for volume rendering the cube and mak- 
ing the PdBI movie. The observations in this paper were 
carried out with the IRAM Plateau de Bure Interferom- 
eter and the Expanded Very Large Array. PdBI is sup- 
ported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). The EVLA is a facility of NRAO, operated 
by AUI, under a cooperative agreement with the NSF. 



REFERENCES 



Andrews, B. H. & Thompson, T. A. 2011, ApJ, 727, 97 
Blain, A. W., Smail, I., Ivison, R. J., & Kneib, J.-P. 1999, 

MNRAS, 302, 632 
Blitz, L. & Rosolowsky, E. 2004, ApJL, 612, L29 
— . 2006, ApJ, 650, 933 

Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F., & Blitz, 

L. 2008, ApJ, 686, 948 
Bonnell, I. A., Clarke, C. J., & Bate, M. R. 2006, MNRAS, 368, 

1296 

Bournaud, F. & Elmegreen, B. G. 2009, ApJL, 694, L158 
Bryant, P. M. & Scoville, N. Z. 1996, ApJ, 457, 678 
Chieze, J. & Pineau Des Forets, G. 1987, A&A, 183, 98 
Chieze, J. P. 1987, A&A, 171, 225 
Courteau, S. 1997, AJ, 114, 2402 

Daddi, E., Bournaud, F., Walter, F., Dannerbauer, H., Carilli, 
C. L., Dickinson, M., Elbaz, D., & Morrison G. et al. 2010, 
ApJ, 713, 686 

Danielson, A. L. R., Swinbank, A. M., Smail, I., Cox, P., Edge, 
A. C, Weiss, A., Harris, A. I., & Baker, , A. J. et al. 2011, 
MNRAS, 410, 1687 

Downes, D. & Solomon, P. M. 1998, ApJ, 507, 615 

Ebeling, H., Edge, A. C, & Henry, J. P. 2001, ApJ, 553, 668 

Elmegreen, B. G. 1989a, ApJ, 338, 178 

— . 1989b, ApJ, 344, 306 

— . 1993, ApJ, 411, 170 

Elmegreen, B. G. & Efremov, Y. N. 1997, ApJ, 480, 235 
Elmegreen, B. G., Elmegreen, D. M., Fernandez, M. X., & 

Lemonias, J. J. 2009, ApJ, 692, 12 
Engel, H., Tacconi, L. J., Davies, R. I., Neri, R., Smail, I., 

Chapman, S. C, Genzel, R., & Cox, P. et al. 2010, ApJ, 724, 

233 

Forster Schreiber, N. M., Genzel, R., Bouche, N., Cresci, G., 

Davies, R., Buschkamp, P., Shapiro, K., & Tacconi, L. J., et al. 

2009, ApJ, 706, 1364 
Gao, Y. & Solomon, P. M. 2004, ApJS, 152, 63 
Genzel, R., Burkert, A., Bouche, N., Cresci, G., Forster 

Schreiber, N. M., Shapley, A., Shapiro, K., & Tacconi, L. J., et 

al. 2008, ApJ, 687, 59 



Genzel, R., Newman, S., Jones, T., Forster Schreiber, N. M., 
Shapiro, K., Genel, S., Lilly, S. J., & Renzini, A. et al. 2011, 
ApJ, 733, 101 

Genzel, R., Tacconi, L. J., Gracia-Carpio, J., Sternberg, A., 

Cooper, M. C, Shapiro, K., Bolatto, A., & Bouche, N. et al.,. 

2010, MNRAS, 407, 2091 
Greve, T. R., Papadopoulos, P. P., Gao, Y., & Radford, S. J. E. 

2009, ApJ, 692, 1432 
Harris, W. E. & Pudritz, R. E. 1994, ApJ, 429, 177 
Heyer, M. H. & Brunt, C. M. 2004, ApJL, 615, L45 
Honma, M., Sofue, Y., & Arimoto, N. 1995, A&A, 304, 1 
Ivison, R. J., Papadopoulos, P. P., Smail, I., Greve, T. R., 

Thomson, A. P., Xilouris, E. M., & Chapman, S. C. 2011, 

MNRAS, 412, 1913 
Ivison, R. J., Swinbank, A. M., Swinyard, B., Smail, I., Pearson, 

C. P., Rigopoulou, D., Polehampton, E., & Baluteau, J.-P.. et 
al. 2010, A&A, 518, L35 

Jappscn, A., Klessen, R. S., Larson, R. B., Li, Y., & Mac Low, M. 

2005, A&A, 435, 611 
Jones, T. A., Swinbank, A. M., Ellis, R. S., Richard, J., & Stark, 

D. P. 2010, MNRAS, 404, 1247 

Julio, E., Kneib, J.-P., Limousin, M., Eh'asdottir, A., Marshall, 

P. J., & Verdugo, T. 2007, New Journal of Physics, 9, 447 
Kassiola, A. & Kovner, I. 1993, ApJ, 417, 450 
Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 
Keto, E., Ho, L. C, & Lo, K. 2005, ApJ, 635, 1062 
Keto, E. R. & Myers, P. C. 1986, ApJ, 304, 466 
Klessen, R. S. 2004, ApSS, 292, 215 

Klessen, R. S., Spaans, M., & Jappsen, A. 2007, MNRAS, 374, 
L29 

Kneib, J.-P., Ellis, R. S., Smail, I., Couch, W. J., & Sharpies, 

R. M. 1996, ApJ, 471, 643 
Kravtsov, A. V. & Gnedin, O. Y. 2005, ApJ, 623, 650 
Krumholz, M. R. & McKee, C. F. 2005, ApJ, 630, 250 
Larson, R. B. 1981, MNRAS, 194, 809 
— . 2005, MNRAS, 359, 211 



Swinbank et al. 



15 



Law, D. R., Steidel, C. C, Erb, D. K., Larkin, J. E., Pettini, M., 

Shapley, A. E., & Wright, S. A. 2009, ApJ, 697, 2057 
Lehnert, M. D., Nesvadba, N. P. H., Tiran, L. L., Matteo, P. D., 

van Driel, W., Douglas, L. S., Chemin, L., & Bournaud, F. 

2009, ApJ, 699, 1660 
Li, Y., Klcssen, R. S., & Mac Low, M. 2003, ApJ, 592, 975 
Mao, R. Q., Henkel, C, Schulz, A., Zielinsky, M., Mauersberger, 

R., Storzer, H., Wilson, T. L., & Gensheimer, P. 2000, A&A, 

358, 433 

Padoan, P. & Nordlund, A. 2002, ApJ, 576, 870 
— . 2011, ApJ, 730, 40 

Papadopoulos, P. P. & Pelupessy, F. I. 2010, ApJ, 717, 1037 
Papadopoulos, P. P. & Seaquist, E. R. 1998, ApJ, 492, 521 
Papadopoulos, P. P., Thi, W.-F., Miniati, F., & Viti, S. 2011, 

MNRAS, 414, 1705 
Plume, R., Jaffe, D. T., Evans, II, N. J., Martin-Pintado, J., & 

Gomez-Gonzalez, J. 1997, ApJ, 476, 730 
Richard, J., Pei, L., Limousin, M., Julio, E., & Kneib, J. P. 2009, 

A&A, 498, 37 

Richard, J., Smith, G. P., Kneib, J. -P., Ellis, R. S., Sanderson, 

A. J. R., Pei, L., Targett, T. A., & Sand, D. J. et al. 2010, 

MNRAS, 404, 325 
Riechers, D. A., Carilli, L. C, Walter, F., Weiss, A., Wagg, J., 

Bertoldi, F., Downes, D., & Henkel, C. et al. 2011a, ApJL, 733, 

Lll 

Riechers, D. A., Hodge, J., Walter, F., Carilli, C. L., & Bertoldi, 

F. 2011b, arXiv/1105.4177 
Rosolowsky, E. & Blitz, L. 2005, ApJ, 623, 826 
Sakamoto, K., Wang, J., Wiedner, M. C, Wang, Z., Peck, A. B., 

Zhang, Q., Petitpas, G. R., & Ho, L., et al. 2008, ApJ, 684, 957 
Schaye, J. & Dalla Vecchia, C. 2008, MNRAS, 383, 1210 
Schweizer, F., Miller, B. W., Whitmore, B. C, & Fall, S. M. 1996, 

AJ, 112, 1839 



Scovillc, N. 2004, in Astronomical Society of the Pacific 

Conference Series, Vol. 320, The Neutral ISM in Starburst 

Galaxies, cd. S. Aalto, S. Huttemeister, & A. Pedlar, 253 — h 
Shapiro, K. L., Genzel, R., & Forster Schreiber, N. M. 2010, 

MNRAS, 403, L36 
Smail, I., Swinbank, A. M., Richard, J., Ebeling, H., Kneib, J. -P., 

Edge, A. C, Stark, D., Ellis, R. S., Dye, S., Smith, G. P., & 

Mullis, C. 2007, ApJL, 654, L33 
Smith, G. P., Ebeling, H., Limousin, M., Kneib, J.-P., Swinbank, 

A. M., Ma, C.-J., Jauzac, M., & Richard, J., et al. 2009, ApJL, 

707, L163 

Solomon, P. M., Downes, D., Radford, S. J. E., & Barrett, J. W. 

1997, ApJ, 478, 144 
Spaans, M. & Silk, J. 2000, ApJ, 538, 115 
Spcrgel, D. N., Verde, L., Peiris, H. V., Komatsu, E., Nolta, 

M. R., Bennett, C. L., Halpern, M., & Hinshaw, G., et al. 2003, 

ApJS, 148, 175 

Swinbank, A. M., Smail, I., Longmore, S., Harris, A. I., Baker, 

A. J., De Breuck, C, Richard, J., & Edge, A. C. et al. 2010, 

Nature, 464, 733 
Tacconi, L. J., Genzel, R., Ncri, R., Cox, P., Cooper, M. C, 

Shapiro, K., Bolatto, A., & Bournaud, F., et al. 2010, Nature, 

463, 781 

Tacconi, L. J., Genzel, R., Smail, I., Neri, R., Chapman, S. C, 
Ivison, R. J., Blain, A., & Cox, P., ct al. 2008, ApJ, 680, 246 
Toomre, A. 1964, ApJ, 139, 1217 

Walter, F., Weiss, A., & Scoville, N. 2002, ApJL, 580, L21 
WeiB, A., Neininger, N., Huttemeister, S., & Klein, U. 2001, 
A&A, 365, 571 

WeiB, A., Walter, F., & Scoville, N. Z. 2005, A&A, 438, 533 
Wilson, C. D., Scoville, N., Madden, S. C, & Charmandaris, V. 

2003, ApJ, 599, 1049 
Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, 

A. G. G. M. 2003, ApJ, 587, 278 



